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Abstract 

We introduce a generalization of the (/-analysis, which provides a novel non-parametric tool 
for the description and detection of log-periodic structures associated with discrete scale invari- 
ance. We use this generalized q analysis to construct a signature called the (H, ^-derivative 
of discrete scale invariance, which we use to detect the log-periodicity in the energy release 
rate and its cumulative preceding the rupture of five pressure tanks made of composite carbon- 
matrix material. We investigate the significance level of the spectral Lomb periodogram of the 
optimal (H, q)-derivative. We confirm and sttengthen previous paramettic results that the en- 
ergy release rate and it cumulative exhibit log-periodicity before rupture. However, our tests to 
use this method as a scheme for the prediction of the critical value of the stress at rupture are 
not encouraging. 

1 Introduction 

The fracture of materials is a catastrophic phenomenon of considerable technological and scientific 
importance. However, a reliable identification of precursory signatures of impending failure is lack- 
ing in most cases. A notable exception has been found [Jlj, ^] in the analysis of acoustic emissions 
recorded during the pressurization of spherical tanks of kevlar or carbon fibers pre-impregnated in a 
resin matrix wrapped up around a thin metallic liner (steel or titanium) fabricated and instrumented 
by Aerospatiale-Matra Inc. (now EADS). A recent thorough analysis ^ of the seven acoustic 
emission recordings of seven pressure tanks, that was brought to rupture, has unambiguously char- 
acterized the acceleration of acoustic energy rate dE/dt and found it to be in agreement with a 
power law "divergence" expected from the critical point theory proposed in Ref. [Q]. In addition, 
strong evidence of log-periodic corrections was found that quantify the intermittent succession 
of accelerating bursts and quiescent phases of the acoustic emissions on the approach to rupture. 
Ref. [^] also proposed an improved model accounting for the cross-over from the non-critical to 
the critical region close to the rupture point exhibits interesting predictive potential. The critical 
rupture concept, confirmed by other experiments opens the road towards industrial applications 
involving heterogeneous materials such as fiber composites, rocks, concrete under compression and 
materials with large distributed residual stresses 

However, a time-to-failure behavior following a power law dE/dt oc (t c — t)~ a does not provide 
a reliable and unique signature: fits of noisy data by such power laws are notoriously unreliable; for 
instance, an error of 1% in the determination of t c usually leads to errors of tens of percent for the 
exponent a. In addition, the determination of t c is very sensitive to the presence of noise. In order 
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to improve the determination of t c , the existence of log -periodic oscillations have been found useful 
^ [7[ |8|] and have been used for the implementation of prediction schemes |2], |3]] with reasonable 
success. 

The demonstration of the existence of log-periodic corrections to the power law is of both 
fundamental and practical interest. From a fundamental point of view, log-periodicity signals a 
spontaneous hierarchical organization of damage with an approximate geometrical set of character- 
istic scales. A possible mechanism involves a succession of ultraviolet instabilities of the Mullins- 
Sekerka type ^ (see also [10] for a review). More generally, the presence of log -periodicity signals 
the partial breaking of continuous scale invariance into discrete scale invariance, which requires that 
the underlying field theory be non-unitary flU|]. From a practical view point, log -periodicity may 
help locking in the fit on experimental data to obtain a better precision on the recovery of the critical 
rupture time t c [j^HUl. 

However, most of the evidence of log-periodicity in rupture results from parametric fits of the 
experimental or numerical data by a log-periodic power law formula, except for Ref. [Q] which 
introduced a "canonical" averaging method to extract the log-periodic signal directly. Parametric 
fits suffer from two problems: (1) the formula cannot avoid some simplification which for instance 
omits the presence of harmonics and/or other structures; (2) parametric fits are delicate due to 
possible degeneracies and, in addition, their statistical significance (i.e., added value) is difficult to 
estimate. It is thus important to develop further non-parametric tests. This is our goal here to present 
a novel non-parametric method, that turns out to be very powerful in identifying log-periodicity in 
noisy data. 

Our method is based on the concept of a ^-derivative, the inverse of the Jackson (/-integral [12], 
which is the natural tool [13, 14] for describing discrete scale invariance. Indeed, (/-derivatives can 
be identified with the generators of fractal and multifractal sets with discrete dilatation symmetries 
[14]. Nowhere differentiable functions which characterize fractal or multifractal sets turn out to 
be perfectly well behaved under the (/-derivative. Discrete renormalization group equations, whose 
general mathematical solutions are power laws with complex exponents (and hence exhibit log- 
periodicity), can be seen as nothing but Jackson (/-integrals of regular functions of the decimated 
degrees of freedom. Jackson (/-integrals constitutes the natural generalization of regular integrals 
for discretely self-similar systems [14]. The way that Jackson (/-integral can be related to the free 
energy of a spin system on a hierarchical lattice was explained in [|13|]. 

In section 2, we introduce the (/-derivative and generalize it to take into account of anomalous 
scaling. We discuss the main properties of the generalized (/-derivative that will be useful for our 
analysis of this paper. Section 3 presents our analysis with the generalized (/-derivative of the 
acoustic emission data, both for the energy release rate and for the cumulative energy release, 
obtained during the pressurization of spherical tanks of kevlar or carbon fibers pre-impregnated 
in a resin matrix wrapped up around a thin metallic liner (steel or titanium). For comparison, we 
use exactly the same set of seven acoustic emission recordings of the seven pressure tanks used in 
the previous study reported in [^]. Our new results confirm strongly the existence of log-periodicity 
with a much enhanced confidence. Section 4 tests a scheme using the generalized (/-derivative for 
prediction purpose. Here, we find disappointing results: the parametric approach in ^ turns out to 
be more powerful. 



2 The generalized (/-analysis and log-periodicity 
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2.1 Definition 



Let us take some q G (0, 1) U (1, oo). The (/-derivative of an arbitrary function f(x) is defined as 

o,m = HplM . (i) 

(1 - q)x 

For q — > 1, the definition ([[]) recovers the usual definition of a derivative. 

For g ^ 1, D q f(x) is more than just a derivative: it compares the relative variations of f(x) 
and of x when x is magnified by the finite factor q. It is thus intuitive that the (/-derivative tests the 
scale invariance property of the function f(x). As we said above, it was actually shown by Erzan 



and Eckmann [ ]14[ ] that the (/-derivative is the natural tool for describing discrete scale invariance, 
since a fixed finite q compares f(x) with f(qx) at x magnified by a fixed factor, and thus it also 
compares f(qx) with f(q 2 x), f(q 2 x) with f(q 3 x), etc. When x is taken as the distance from a 
critical point, D q f(x) thus quantifies the discrete self-similarity of the function f(x) in the vicinity 
of the critical point. From the definition ([[]), it is clear that 

D 1/q f(x) = D q f(x/q) . (2) 

It is thus enough to study D q f(x) for q G (0, 1) to derive its values for all q's. 

The necessary and sufficient condition for a function f(x) of order ip to be homogeneous is 

q — 1 x 

This expression suggests the introduction of a generalized (/-derivative that we call (H, q) -derivative, 
such that the dependence in x of D q f(x) disappears for homogeneous functions, for the choice 
H = ip. Consider therefore the following definition 

such that D q =1 f(x) recovers the standard (/-derivative D q f(x). For a power law function f(x) = 
Bx m , D q H=m [Bx m ] = B(l — q m ) /(l — q) m is constant. For a statistically homogeneous function 

f{x) = Bx m , D^ =1 f(x) = constant. 

The generalized (H, q) -derivative has two control parameters: the discrete scale factor q devised 
to characterize the log-periodic structure and the exponent H introduced to account for a possible 
power law dependence, i.e., to trends in log-log plots. 

2.2 Application to log-periodic functions 

Since Erzan and Eckmann [ |14| ] showed that the (/-derivative is the natural tool for describing discrete 
scale invariance, it is natural to study the properties of the (H, q) -derivative of the simplest func- 
tion exhibiting discrete scale invariance, namely a power law decorated by a simple log-periodic 
function: 

y(x)=A- BT m x m + Cr m x m cos (u> In x) , (5) 

where the presence of a phase <\> = lu In r has been absorbed in the definition of x, < C < B and 
to = 2itf. This equation, where x is interpreted as the normalized distance x = (p c — p)/r to a 
critical point p c , has been used in several works to describe material rupture [jl|, 0, |L |h, precursory 
patterns of large earthquake s [|T3[ ], rock bursts [pJ], aftershocks [17, 18], and speculative bubbles 
preceding financial crashes [|19|]. 
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The (H, q) -derivative of y(x) is 

D«y(x)=x m - H [B , + C'g(x)}, (6) 

where 

Br m (l-q m ) Cr m 

and 

= Ci cos(a; In x) + C2 sin(ct> In x) , (8) 

with C\ = 1 — q rn cos(w In g) > and C2 = sin(a; In g). The special choice H = m gets rid of the 
power law and the (H, q) -derivative D^y{x) becomes the pure log-periodic function g(x). 

For the special choices q = e r n2 ' K l lJJ where n is a positive integer corresponding to choosing q 
equal to one of the preferred scaling factor of the log-periodic function, we obtain C\ = 1 — q m and 
C2 = 0. Therefore, the (H, g) -derivative offers a novel approach for detecting the preferred scaling 
factors of the discretely scale invariant function by a measure of its phase: those values of q such 
that the phases of the (H = m, g)-derivative is that of a pure cosine should qualify as the preferred 
scaling factors. Such phases can be for instance measured by the Hilbert transform. Here, we do 
not persue this possibility which wll be explored in another presentation. 

In the Appendix A, it is shown that g(x) is extremal at x rn solution of 

uln(x m ) = nir + arctan(C2/Ci) , (9) 
where n = 0, ±1, ±2, • • •, and that the extreme values of g(x) are 

g(x m ) = ±Jc[+Ci . (10) 

The amplitude of D^y(x) is then 

A = xZ~ H C yf(% + C| , (11) 
while the successive extreme values are 

D m = (B' ± C'y/Cf + Cf) (x m ) m - H , (12) 

where B' and C are defined in (^). By fixing H close to m, one can in principle obtain the 
amplitude A m to be constant as a function of x. This value H = m should provide theoretically the 
most significant log-periodic component, quantified for instance by the largest peak of the Lomb 
periodogram. However, in practice, the noise embedded in the data may distort the log-periodic 
oscillations and the most significant log-periodicity may occur for H / m. The introduction of H 
affords a convenient detrending scheme. 

Figure [l] shows y{x) defined in Eq. (||) as a function of the distance x to the critical point with 
A = 1260, B = 300, C = 6, m = 0.3, u> = 5.4 and <j) = 0. Figure |2] shows its generalized 
g-derivative with q = 0.5 for H = 0.2, H = 0.3 and H = 0.4, respectively. This generalized 
g-derivative has been calculated by using the incorrect assumption that the critical point is at x = 1 
in order to also show the distortion resulting from an error in the determination of the critical point. 
This distortion becomes important when x is not large compared to 1. We observe that the amplitude 
of the oscillations of D^y(x) increases with H when going towards the critical point x = 0, in 
agreement with the prediction (p4|). 

In the next section, using the (H, g)-analysis, we test for the presence of log-periodicity in the 
energy release rate and its cumulative obtained in the experimental recordings of the seven pressure 
tanks used in the previous study reported in [§]. 
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3 Investigation of the significance of log-periodicity in acoustic emis- 
sion energy release using the generalized (if, q) -analysis 

We use the same notations to label the data sets as used in Ref. [||]. In this analysis, we use the 
true value p c of the pressure at which rupture occurred in order to define the pressure-to-failure 
p c — p quantifying the distance to the critical rupture point. We follow ^ and use the acoustic 
emission data from the pressure at which a noticeable acceleration in the cumulative energy release 
takes place. This choice is not crucial at all, since the (H, g)-analysis is not sensitive to the points 
far from the critical point. We also exclude the last six points nearest to p c because they contain 
the largest noise and may suffer from finite size effects that may lead to serious distortions. In 
our analysis, we exclude data sets of pressure tanks # 5 and # 7, because for these two data sets, 
Piast <C p c (the last data point is for a pressure far below the critical rupture value). This leads to 
few oscillations and low statistical significance. We thus are left to analyze five data sets # 1, # 2, 
#3, #4 and #6. 



3.1 Energy release rate 

In each analysis, q ranges from 0.1 to 0.9 with spacing 0.1, while H takes values from —0.9 to 
0.9 with spacing 0.1. This defines 9 x 19 = 171 parameter pairs (H, q). For each parameter pair 
(H,q), we calculate the generalized (/-derivative of the energy release rate. We then perform a 
spectral analysis of the generalized (/-derivative using the Lomb periodogram method [^0|], in order 
to test for the statistical significance of possible log-periodic oscillations. For each (H, q) pair, 
the highest peak Pn(H, q) and its associated logarithmic angular frequency to(H, q) in the Lomb 
periodogram are obtained. The first criterion used to identify a log-periodic signal is the strength 
of the Lomb periodogram analysis, i.e., the height of the spectral peaks. 



The results for the five data sets # 1, # 2, # 3, # 4 and # 6 are shown in Figs. [| to [22| Figures 
|], H [l2|, [if] and 2C show the dependence of the highest peak P/v (H, q) in each Lomb periodogram as 
a function of H and q. Figs. |3],[7, ll|, 15 and [l^ give the associated logarithmic angular frequencies 
u>(H, q). The other figures show the generalized (/-derivative for the pair (H, q) giving the highest 
Lomb peak and the associated Lomb periodogram. 

In order to interpret the results presented in the figures ||]to 22, we need to explain the criteria 



that we have used. First, to make the description simpler and more geometrical, as shown in Fig. |[ 
we first classify all pairs of (H, q) into three classes: W (wedge or wall), P (platform) and B 
(bottom of valley or basin). 

In Fig. |], the pairs (H, q) near (0.1, 0.1) gives the largest Pat, which implies the most significant 
log-periodic oscillations in the generalized q derivative. However, the associated u> is about 3 ~ 4 
as shown in the region B of Fig. ||. The problem is that this value is dangerously close to the most 
probable (log-periodic) angular-frequency u> mp = 3.6, resulting solely from the most probable 
noise. It was indeed shown [ |l7| ] that noise decorating power laws may lead to artifactual log- 
periodicity with a most probable frequency corresponding roughly to 1.5 oscillations over the whole 
range of analysis. In the present context, this defines the most probable (log-periodic) angular 
frequency oj mp by the following formula 

u m P = 2 * X L5 (13) 

ln(p c - p min ) - ln(p c - p max ) ' 

where the acoustic emission signal is recorded from the pressure p m j n to p max - In the case of 
experiment # 1, this leads to uj mp = 3.6. Thus, a value of uo in the range 3 ~ 4 corresponds to 
approximately 1.5 oscillations in the plot of the generalized (/-derivative as a function of ln(p c — p) 
for the whole range of pressure. This is not sufficient. We should thus exclude these pairs of 
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(H,q). This defines our second criterion for qualifying the existence of log -periodicity: pairs 
of (H, q) with oj(H, q) < uj mp should be discarded, because there is a non-negligible probability 
that the observed log-periodicity may result from noisy fluctuations around the power law and is 
thus spurious. Ideally, we should estimate the probability that random noise, of several plausible 
standard distributions, creates a false alarm that a periodicity (or log-periodicity) is found in the 
(H, q) -derivative of the signal. This has been done in a systematic way for a large variety of noise, 
without and with long-range correlation [21]. However, it is difficult to identify what should be 
the correct null hypothesis of the noise decorating the generalized ^-derivative. For the sake of 
simplicity, we thus stick to the simple rule called second criterion just discussed. 

A log-periodic angular frequency which is too small may result from noise, as we just said. 
Similarly, a too high log-periodic angular frequency is also the signature of noise, simply because 
it is easy to fit noisy data with many oscillations. We should thus also discard the pairs of (H, q) 
whose lo are too large. Ref. [J3j] used this criterion to discard solutions with uj > 14. This gives us a 
third criterion. It is not always obvious to fix the value of the threshold angular frequency beyond 
which solutions are rejected as noise. Fortunately, this third criterion can often be obtained to 
follow from the first criterion, since the Lomb peak corresponding to a large angular frequency 
is usually low [21]. We apply this criterion and hence exclude (H, q) pairs in the regions of Wi, 
W2 and W3 shown in Fig. H. 

A physically meaningful optimal pair of (H, q) is thus obtained. For data set # 1, it corresponds 
to (H, q) = (—0.9, 0.5), which is indicated by arrows in Fig. || and Fig. |j. It lies within a platform 
shown in Fig. [| We also show in Fig. || the generalized (/-derivative and in Fig. |6] its corresponding 
Lomb periodogram. The log-periodic oscillations are found to be very significant. Under the null 
hypothesis of independent Gaussian noise decorating the signal, the false-alarm probability of the 
log-periodic component is k 0.06%. The false-alarm probability for the null hypothesis of inde- 
pendent heavy-tailed noises (say, Levy stable noise and power-law noise) and weakly correlated 
noises (say, Garch(l,l) noise and fractional Gaussian noise (fGn) with the Hurst index less than 
0.5) is even lower [p|]. If we assume a strongly correlated noise such as fractional Gaussian noise 
with the Hurst index greater than 0.5, the false-alarm probability increases: for instance, a fractional 
Gaussian noise with Hurst index of 0.62 leads to a false-alarm probability of 1% [21]. 

As a fourth criterion, if the majority of the pairs (H, q) have similar oscillatory behavior (i.e., 
similar uj), we take that this indicates a robust log-periodic signal and that the optimal pair should 
be within these pairs with high probability. This excludes Wi, W2, W3 and W4 in Fig. || and 
strengthens the choice of the previous criteria. This fourth criterion is sometimes ambiguous to 
apply in practice. However, all experimental systems have their optimal pair (H, q) in well-defined 
platforms P, except for experiment # 3 which presents its optimal (H, q) in Wi. 

Table @ summurizes the results of the (H, g)-analysis on the energy release rate of data sets 
# 1, # 2, # 3, # 4 and # 6. For each data set, the generalized (/-analysis was performed between 
(Pmin)Pmax)> with a number of points varying between 55 and 164 as indicated in the column 
"Points." p c is the true critical pressure at rupture. Pjy is the optimal value of the Lomb peak 
height, and lo is the corresponding logarithmic angular frequency. The column "Gaussian" presents 
the false-alarm probability of the Lomb peak under the null hypothesis of independent Gaussian 
noise. The column "fGn" evaluates the value of the Hurst index of a fractional Brownian noise 
which would give a false-alarm probability 1% to get the same peak as in our analysis. 

These results give a very strong confidence that there is a genuine log-periodicity of the energy 
release rate before rupture in the data sets # 1, # 2, # 3, # 4 and # 6. This confirms and strengthens 
Ref. [g], which claimed that a pure power law fails to parameterize the data but a log-periodic 
formulae does a good job. 
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3E-13 
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9E-3 


0.52 


#6 


164 


451.5 


723.5 
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(0.1,0.2) 


15.8 


6.2 
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2.9 


3.2 


2E-5 


0.71 



Table 1: Summary of the results of the (H, q) -analysis on the energy release rate of data sets 
# 1, # 2, # 3, # 4 and # 6. For each data set, the generalized g-analysis was performed between 
(PmimPmax) with a number of points given in the column "Points." p c is the true critical pressure 
of rupture. The column (H, q) lists the optimal pairs. Pjv is the optimal value of the Lomb peak 
height. u> and u 1 are the corresponding logarithmic angular frequency and the angular frequency 
corresponding to the second highest peak. N osc is the number of oscillations. The column "Gaus- 
sian" presents the false-alarm probability of the Lomb peak under null hypothesis of independent 
Gaussian noise. The column "fGn" gives the Hurst index of a fractional Gaussian noise that would 
give a false-alarm probability 1% to obtain the same Lomb peak as in the signal. 



3.2 Cumulative energy release 



We follow the same procedure as in the previous section to analyze the cumulative energy release 
of the same five experimental data sets discussed previously. We use slightly different truncations 
in this analysis of the cumulative energy release. The results obtained for the five data sets are 
shown in Figs. ^ to 42. Figures 24, 28, [32|, 36 and 40 show the dependence of the highest peak 

ill, |35| and M are 



Pn(H, q) in each Lomb periodogram as a function of H and q. Figs. 23, 27 
plots of the associated logarithmic angular frequencies uj{H. 



The generalized (/-derivative of 
the cumulative energy release corresponding to the optimal pairs of (H, q) for the five data sets 
are illustrated respectively in Figs. 25, 29, 33, 37 and |1] with their Lomb periodograms shown in 
Figs. ^[^,[M[[|8| and ||. 



In Fig. 24, the largest P/v's are in the region B. However, the corresponding frequencies are 



very low as shown in Fig. [23| which gives the most probable frequency. Hence we excludes region 
B in accordance with the second criterion. The optimal pair in the platform P is (H = — 0.5, q = 
0.6). One can observe a platform of u> in Fig. ^3], which is consistent with the fourth criterion. In 
Fig. the log-periodic oscillations of the (H, g)-derivative are clearly visible, and correspond to 
the high Lomb peak shown in Fig. The log-periodic oscillations become even strongly and more 
significant if more points are discarded at the low and high pressure ends of the data set. 

Data set # 2 is more noisy and a major part of the pairs (H, q) give very low frequencies close 
to uj mp , as shown in Fig. ^7|, which lead us to exclude all the local maxima shown in Fig. 28. The 
optimal pair (0.3,0.1) is located at the wedge (W3) indicated by an arrow (in contradiction the 
fourth criterion). The log-periodicity is quite apparent in Fig. |9[ Correspondingly, its Lomb peak 
shown in Fig. [3(] is very significant. The data sets # 3 and # 4 are similar to # 2. For # 6, the 
optimal pair (—0.3, 0.8) located on the platform shown in Fig. ^9] has a local maximum shown in 
Fig. ^ We summarize all these results in Table |2[ 



3.3 Comparison between energy release rates and their cumulative 

Comparing Table [T] and Table ^ we find that the cumulative energy releases present a more signif- 
icant log-periodic component. Taking the cumulative corresponds to perform a low-pass filter that 
smooths out significantly the noise and usually provides better signals with higher signal-to-noise 
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17.5 
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0.79 


#6 
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(-0.3,0.8) 


15.7 


12.7 


5.4 


3.9 


4.9 


1E-5 


0.71 



Table 2: Summary of the results of the (H, g)-analysis on the cumulative energy release of data sets 
# 1, # 2, # 3, # 4 and # 6. The meaning of the columns are the same as in Table |. 



ratio. Of course, taking the cumulative also reduces the amplitude of the underlying log-periodic 
oscillations. There is thus a compromise. As an illustration, let us compare Fig. 13 and Fig. |33|. It is 
clear that the (H, q) -derivative of the energy release rate in Fig. 13 is much more noisy than that of 
the cumulative energy release shown in Fig. [33|. On the other hand, the high-frequency oscillations 
which are clearly visible in Fig. [53| are suppressed in Fig. 33 and can not be captured by the Lomb 
analysis. Thus, the log-frequency with largest peak power of the energy rate (/ = 1.97) is much 
higher than that of the cumulative energy (/ = 0.72). We note however that the second highest 
peak in Fig. 33 corresponds to the highest peak of the energy release rate. 

In general, we find that the cumulative energy release keeps the same underlying log-periodic 
structure as found in the energy release rate. We observe that data sets # 1, #2 and #4 have 
very similar logarithmic frequencies for the energy release rate and their cumulative. Most Lomb 
periodograms exhibit harmonics of each other. This close correspondence obtained in the (H, q)- 
analysis provides a better characterization of the log-periodic nature in the system than obtained 
previously with the parameterization approach of fitting the data sets with a log-periodic function 

§■ 

Finally, we note that all data sets exhibit a fundamental value of the logarithmic angular fre- 
quency which is compatible with the value to = 5.8 ± 0.9. In some cases, it is its second harmonic 
that gives the strongest peak of the Lomb periodograms. This value corresponds to a preferred 
scaling ratio A = e 27r/uj = 3.0 ± 0.5. 



4 Post-diction of ruptures 

The (H, g)-analysis has shown its power for detecting log-periodicity, conditioned on our knowl- 
edge of the critical pressure p c at rupture. It is natural to ask whether it can be extended to provide 
advanced prediction of p c . For carrying these tests, we use the cumulative energy release which pro- 
vided the strongest log-periodic signal in the analysis reported in previous sections. Our strategy 
is to use each data up to a maximum pressure p max < Pc, assume some value for p c and perform 
the same analysis as in section 3. i.e., for a given presumed critical points p c , we determine the 
optimal pair (H,q). For each presumed p c , we thus obtain the Lomb peak height Pjy(H,q) and 
its associated logarithmic angular frequency oj(H, q) as a function of H and q. In order to find the 
optimal (H, q) corresponding to the maximum Pn{H, q), we add the second criterion to exclude 
those (H, q) with u>(H, q) < u> mp . Having determined the optimal Pn{p c ) and u(p c ) as functions 
of the presumed critical pressure p c , we determine the predicted critical pressure p c by the condition 

P N {p c ) = max{P N (p c )}, (14) 

Pc 

without further constraints on u)(p c ). 
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Points 


Pmin 


Pmax 


to 


Pc 


Pc Pmax 


Pc l) 


.2 
Pc 


#1 


139 


453.5 


705.5 


5.5 


713 


1.1% 


722.5(734) 


722.5(761) 


#2 


149 


452.5 


662.5 


5 


673 


1.6% 


680.5(681) 


686.5(700) 


#3 


65 


671.5 


759.5 


4.5 


764 


0.6% 


763.5(765) 


767.5(772) 


#4 


56 


671.5 


746.5 


5.5 


756 


1.3% 


755.5(755) 


755.5(763) 


#6 


74 


614.5 


726.5 


6.4 


734 


1.0% 


742.5(744) 


758.5(759) 



Table 3: Post-diction using the (H, g)-analysis on the cumulative energy release of tanks # 1, # 2, 
# 3, # 4 and # 6. Columns pi 1 ^ and p£ 2 ^ list the predicted critical pressure p c and its corresponding 
upper threshold in the parenthesis with the constraint ([16]) for k = 4/3 and k = 1 respectively. The 
meaning of the other columns are the same as in Table |]. 



In practice, we analyze the function P^{p c ;H,q) of the three variables p c , H and q in the 
following sequence: we fix a value for p c and explore the plane (H,q). We then change p c and redo 



the exploration of the plane (H, q), etc. Fig. |3| through 47 give the optimal Lomb peak height Pjy 
as a function of the presumed critical pressure p c . It is obvious that this scheme does not give a good 
predictive skill, as the most important pattern observed in these figures is the systematic increasing 
trend, tending to reject the predicted p c towards too large values. 

We thus attempt to improve this prediction skill by requesting that the predicted p c should not 
be too far from the last point, i.e., it is nonsensical to hope to predict too far from the "present." To 
implement this idea, we impose u> > kui mp , namely, 

> 1-5/c ^ 

2ir ln(p c - p min ) - ln(p c - p max ) ' 

where k > 1 is a "safety factor." This constraint translates into the following condition for p c : 

. Pmax — Pmin . . 

Pc<PmaxH 3^ • ( 16 ) 

e~ - 1 



This constraint (16) means that there exists an upper threshold beyond which we can not make 
a physically meaningful prediction. Since the left-hand-side of ( |l6| ) is monotonously increasing 
with lo, it is possible to make a prediction much earlier before the critical point, the larger lo is. 
This is natural since large u implies more log-periodic oscillations and thus a stronger log-periodic 
signal. To implement this condition in practice, we take lo = 5.8 as the central value of the log- 
periodic angular frequency. According to this constraint (|T6|), ruptures of tanks # 5 and # 7 are 
unpredictable, since the last point pi as t is too far from the true critical pressure. 

The results are given in Table ||. The columns p^ and p^ list the predicted critical pressure p c 
and its corresponding upper threshold in the parenthesis with the constraint ( |I6| ) for the two choices 
k = 4/3 and k = 1 respectively. The prediction for experiment # 1 is very good. For the other 
experiments, p c is found to be very close to the upper threshold. The predictions for experiments 
# 2, # 3 and # 4 are reasonable while the prediction for case # 6 fails completely. 

Acknowledgments: We are grateful to A. Erzan for a discussion on Jackson's integral and for 
supplying the corresponding references. This work was supported by NSF-DMR99-71475 and the 
James S. Mc Donnell Foundation 21st century scientist award/studying complex system. 
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APPENDIX A 

From the definition (||), the extremal condition dg(x)/dx = gives the extreme points of g{x) 
solutions of 

Ci sin(a; lnz) — C2 cos(o; lnx) = , (17) 

that is, 

to ln(x m ) = nir + arctan(C2/Ci), (18) 
where n = 0, ±1, ±2, • • •. It follows that the extreme values of g(x) are 

g(x m ) = ±y/(% + Cf . (19) 



When ijj ln(x m ) = 2n7r+arctan(C2/Ci), g(x m ) = yCf + C| is the maximum; whencjln(x m ) 



(2n + l)7r + arctan(C2/Ci), g(x m ) = —\jc\ + Cf is the minimum. This is proved by the fol- 
lowing exhaustive classification. 

1. e (2*-i)"-/w < ? < e 2*»r/« with k = 0, ±1, ±2, • • -. 

In this case, C2 < 0. Thus — tt/2 < arctan(C2/Ci) < 0. When u>ln(x m ) = 2nir + 
arctan(C2/Ci), cos^lnxm) > and siri(uj\nx m ) < 0, and thus g{x m ) = \jc\ + Cf 
is the maximum. When wln(x m ) = (2n + l)ir + arctan(C2/Ci), cos(wlna; m ) < and 
sin(o; lnx m ) > 0, and thus g(x m ) = — J Cf + Cf is the minimum. 



2. q = e k7r / w , with k = 0, ±1, ±2, • • •. 

In this case, C2 = 0. Thus arctan(C2/Ci) = and g(x) = C\ cos(o; lnx). When 
u;ln(z m ) = 2nir, g(x m ) = C\ is the maximum. When uj\n(x m ) = (2n + l)7r, g(x m ) = 
— Ci is the minimum. 

3 _ e 2fc7r/ W < ^ < e (2fc+l) 7 r/a; ) w j th fc = Q) ±lj ±2 , • • •. 

In this case, C2 > 0. Thus < arctan(C2/Ci) < 7r/2. When u>ln(x m ) = 2nir + 
arctan(C2/Ci), cos(wlna; m ) > and sin(o;lna; m ) > 0, and thus g{x m ) = \jc\ + C\ 
is the maximum. When uj\n(x m ) = (2n + l)ir + arctan(C2/Ci), cos(wlnx m ) < and 
sin(tj lnx m ) < 0, and thus g(x m ) = — \jc\ + Cf is the minimum. 

The amplitude of D^y(x) is then 

A = C'y/Cf + Cf (x m ) m - H , (20) 
while the local extreme values are 

D m = (b' ± CyjCl + Cf) (x m ) m - H . (21) 
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Figure 1: Plot of y(x) defined by Eq. (Q) as a function of the pressure-to-rupture x with A = 1260, 
B = 300, C = 6, m = 0.3, u = 5.4 and = 0. We generated 120 evenly spaced data points with x 
between 1 and 120, to mimic the cumulative energy release of a real acoustic emission experiment. 




Figure 2: The generalized g-derivative of y(x) shown in Fig. [I] with q = 0.5 for H = 0.2, H = 0.3 
and H = 0.4, respectively. The calculation assumes a critical point at x c = 1 while the synthetic 
function has its genuine critical point at x c = 0. This incorrect value of the critical point in the 
calculation of the generalized (/-derivative is responsible for the distortions observed for x < 1. 
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Figure 3: Dependence of the logarithmic angular frequency co(H, q) of the most significant peak 
in each Lomb periodogram of the (H, q) -derivative of the energy release rate before the rupture of 
tank # 1. The regions Wi, Wj, W3, W4 and B are excluded and the optimal pair (—0.9, 0.5) is 
indicated by an arrow in the platform Pi. 
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Figure 4: Dependence of the height P^{H,q) of the most significant peak in each Lomb peri- 
odogram of the (H, g)-derivative of the energy release rate before the rupture of tank # 1. The 
optimal pair (—0.9, 0.5) is indicated by an arrow in the platform Pi. 
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Figure 5: (H, q) -derivative of the energy release rate before the rupture of tank # 1 as a function of 
the pressure-to-rupture p c — p with q = 0.5 and H = —0.9. 
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Figure 7: Dependence of the logarithmic angular frequency uj(H, q) of the most significant peak 
in each Lomb periodogram of the (H, q) -derivative of the energy release rate before the rupture of 
tank # 2. The wall W is excluded by both the third and the fourth criteria, while B is excluded by 
the second criterion. 




Figure 8: Dependence of the height P^{H,q) of the most significant peak in each Lomb peri- 
odogram of the (H, q) -derivative of the energy release rate before the rupture of tank # 2. 
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Figure 9: (H, q) -derivative of the energy release rate before the rupture of tank # 2 as a function of 
the pressure-to-rupture p c — p with q = 0.1 and H = 0.1. 
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Figure 10: Lomb power of the (H, q) -derivative shown in Fig. |9| 
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Figure 11: Dependence of the logarithmic angular frequency uj(H, q) of the most significant peak 
in each Lomb periodogram of the (H, q) -derivative of the energy release rate before the rupture of 
tank # 3. Regions B and W2, W3, W4 and P are excluded by the second and the third criteria, 
respectively. The fourth criterion is not used in this case and the optimal pair (—0.2, 0.1) is at a wall 
Wi as indicated by an arrow. 
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Figure 12: Dependence of the height P^{H,q) of the most significant peak in each Lomb peri- 
odogram of the (H, q) -derivative of the energy release rate before the rupture of tank # 3. 
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Figure 13: (H, q) -derivative of the energy release rate before the rupture of tank # 3 as a function 
of the pressure-to-rupture p c — p with q = 0.1 and H = —0.2. 
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Figure 14: Lomb power of the (H, g) -derivative shown in Fig. p 
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Figure 15: Dependence of the logarithmic angular frequency u>(H, q) of the most significant peak 
in each Lomb periodogram of the (H, q) -derivative of the energy release rate before the rupture of 
tank # 4. The optimal pair (0.3, 0.4) is in the platform P 2 indicated by an arrow. 
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Figure 16: Dependence of the height Pn(H, q) of the most significant peak in each Lomb peri- 
odogram of the (H, q) -derivative of the energy release rate before the rupture of tank # 4. 
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Figure 17: (i7, g) -derivative of the energy release rate before the rupture of tank # 4 as a function 
of the pressure-to-rupture p c — p with q = 0.4 and H = 0.3. 
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Figure 18: Lomb power of the (H, q) -derivative shown in Fig. 17, The fundamental log-frequency 
/ = 0.76 and its harmonic / = 2.71 are indicated by arrows. 
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Figure 19: Dependence of the logarithmic angular frequency u>(H, q) of the most significant peak 
in each Lomb periodogram of the (H, q) -derivative of the energy release rate before the rupture of 
tank # 6. The optimal pair (0.1, 0.2) is indicated by an arrow. 
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Figure 20: Dependence of the height P^(H, q) of the most significant peak in each Lomb peri- 
odogram of the (H, q) -derivative of the energy release rate before the rupture of tank # 6. 
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Figure 21: (i?, g) -derivative of the energy release rate before the rupture of tank # 6 with respect 
to the pressure-to-rupture p c — p with q = 0.2 and i/ = 0.1. 
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Figure 23: Dependence of the logarithmic angular frequency u(H, q) of the most significant peak in 
each Lomb periodogram of the (H, q) -derivative of the cumulative energy release before the rupture 
of tank # 1. The wedge W and the bottom B are excluded by the second and fourth criteria. The 
optimal pair (—0.5, 0.6) is indicated by an arrow in the platform P. 




Figure 24: Dependence of the height Pjq{H,q) of the most significant peak in each Lomb peri- 
odogram of the (H, g) -derivative of the cumulative energy release before the rupture of tank # 1. 
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Figure 25: (H, q) -derivative of the cumulative energy release before the rupture of tank #1 as a 
function of the pressure-to-rupture p c — p with q = 0.6 and H = —0.5. 
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Figure 26: Lomb power of the (H, g) -derivative shown in Fig. 
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Figure 27: Dependence of the logarithmic angular frequency uj(H, q) of the most significant peak 
in each Lomb periodogram of the (H, q) -derivative of the cumulative energy release before the 
rupture of tank # 2. 
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Figure 28: Dependence of the height P^{H,q) of the most significant peak in each Lomb peri- 
odogram of the (H, g) -derivative of the cumulative energy release before the rupture of tank # 2. 
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Figure 29: (H, q) -derivative of the cumulative energy release before the rupture of tank #2 as a 
function of the pressure-to-rupture p c — p with q = 0.2 and H = 0.3. 
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Figure 30: Lomb power of the (H, g)-derivative shown in Fig. 29. The fundamental log-frequency 
/ = 0.84 and its harmonic / = 1.60 are indicated by arrows. 
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Figure 31 : Dependence of the logarithmic angular frequency u(H, q) of the most significant peak in 
each Lomb periodogram of the (H, q) -derivative of the cumulative energy release before the rupture 
of tank # 3. There are two wedges at q ~ 0.1 and 0.2 excluded according to the fourth criterion, 
around which is the bottom of a basin. The platform with u> > 15 is excluded according to third 
criterion. The optimal pair (—0.2, 0.4) is located within the second platform [—0.2, 0.9] x [0.4, 0.7]. 




Figure 32: Dependence of the height P^(H, q) of the most significant peak in each Lomb peri- 
odogram of the (H, g) -derivative of the cumulative energy release before the rupture of tank # 3. 
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Figure 33: (H, q) -derivative of the cumulative energy release before the rupture of tank # 3 as a 
function of the pressure-to-rupture p c — p with q = 0.4 and H = —0.2. 
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Figure 34: Lomb power of the (H, q) -derivative shown in Fig. [33|. The logarithmic frequencies of 
the three highest peaks are 0.72, 1.57 and 2.61, respectively. 
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Figure 35: Dependence of the logarithmic angular frequency uj(H, q) of the most significant peak 
in each Lomb periodogram of the (H, q) -derivative of the cumulative energy release before the 
rupture of tank # 4. There are two wedges at q = 0.3 excluded according to the fourth criterion. 
The "wall" with u > 15 is excluded according to third criterion. The optimal pair (0, 0.6) is located 
within the platform indicated by an arrow. 




Figure 36: Dependence of the height Pj^(H,q) of the most significant peak in each Lomb peri- 
odogram of the (H, -derivative of the cumulative energy release before the rupture of tank #4. 
The optimal pair (0, 0.6) is indicated by an arrow. 
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Figure 37: (H, q) -derivative of the cumulative energy release before the rupture of tank #4 with 
respect to the pressure-to-rupture p c — p with q = 0.6 and H = 0. 
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Figure 38: Lomb power of the (H, q) -derivative shown in Fig. [37|. The logarithmic frequencies of 
the two highest peaks are 0.87 and 2.79, respectively. 
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Figure 39: Dependence of the logarithmic angular frequency uj(H, q) of the most significant peak 
in each Lomb periodogram of the (H, q) -derivative of the cumulative energy release before the 
rupture of tank # 6. The optimal pair (—0.3, 0.8) is located within a platform. 
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Figure 40: Dependence of the height Pjv (H, q) of the most significant peak in each Lomb peri- 
odogram of the (H, q) -derivative of the cumulative energy release before the rupture of tank # 6. 
The optimal pair (—0.3, 0.8) is indicated by an arrow. 
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Figure 41: (H, q) -derivative of the cumulative energy release before the rupture of tank #6 as a 
function of the pressure-to-rupture p c — p with q = 0.8 and H = —0.3. 
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Figure 42: Lomb power of the (H, q) -derivative shown in Fig. 41. The fundamental log-frequency 
/ = 0.86 and its harmonic / = 2.01 are indicated by arrows. 
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Figure 43: Post-diction of the critical pressure of the rupture of tank # 1. 
indicate the upper threshold of predictable p c estimated by Eq. (16) with k = 
coarse arrows indicate the predicted critical pressures p)t\ p^ and p^\ 



The two fine arrows 
: 4/3 and k = 1. The 
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Figure 44: Post-diction of the critical pressure of the rupture of tank # 2. Same as in figure 
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Figure 45: Post-diction of the critical pressure of the rupture of tank # 3. Same as in figure [13|. 





36 



